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An analytical expression of probability density function (PDF) of accelerations in 
turbulence is derived with the help of the statistics based on generalized entropy 
(the Tsallis entropy or the Renyi entropy). It is revealed that the derived PDF 
explains the one obtained by Bodenschatz et al. in the measurement of fluid 
particle accelerations in fully developed turbulence at R\ = 970. 
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The multifractal analysis of turbulence by the statistics based on the generalized 
entropy of Renyi's or of Tsallis' has been developed by the present authors [ [TJ 
[21HHII1E1E1IZ1E1E]- The Renyi entropy [ [TU] has the extensive character as the 
usual thermodynamical entropy does, whereas the Tsallis entropy [ [TTJ H21 CH] is 
non-extensive. The multifractal analysis belongs to the line of study based on a 
kind of ensemble theoretical approaches that, starting from the log-normal model [ 
H3J E3 EH] , continues with the /3-model [Ej, the p- model [1131111, the 3D binomial 
Cantor set model [ EH] and so on. After a rather preliminary investigation of the 
p-model [HJ, we developed further to derive the analytical expression for the scaling 
exponents of velocity structure function [ElEHHE], and to determine the probability 
density function (PDF) of velocity fluctuations [El El El El and of velocity derivative 
[E] by a self-consistent statistical mechanical approach. 

In this paper, we will derive the formula for the PDF of the accelerations of a 
fluid particle in fully developed turbulence by means of the multifractal analysis. 
With the theoretical PDF, we will analyze the PDF of accelerations at R\ = 970 
(the Taylor microscale Reynolds number) obtained in the Lagrangian measurement 
of particle accelerations that was realized by Bodenschatz and co-workers [ El 122] 
by raising dramatically the spatial and temporal measurement resolutions with the 



T. Arimitsu and N. Arimitsu 



help of the silicon strip detectors. 

We assume that the turbulent flow, satisfying the Navier-Stokes equation 

du/dt + (u-V)u = -V (p/p) + uW 2 u (1) 

of an incompressible fluid, consists of a cascade of eddies with different sizes £ n = S n £o 
where 5 n = 2~ n (n = 0, 1, 2, • • •). The quantities p, p and v represent, respectively, 
the mass density, the pressure and the kinematic viscosity. The acceleration a of a 
fluid particle is given by the substantive time derivative of the velocity: a = du/dt + 
(u ■ V)w. At each step of the cascade, say at the nth step, eddies break up into two 
pieces producing an energy cascade with the energy-transfer rate e n that represents 
the rate of transfer of energy per unit mass from eddies of size £ n to those of size i n +\ 
(the energy cascade model). The Reynolds number Re of the system is given by Re = 
5uq(.q/v = (io/v) 4 ^ 3 with the Kolmogorov scale [I2BJ V = (V 3 /e) 1//4 where e (= eo) 
is the energy input rate to the largest eddies with size £$} Introducing the pressure 
(divided by the mass density) difference 5p n at two points separated by the distance 
£ n , i.e., Sp n = |p/p(«+£ n )—p/p(»)|, and the acceleration a n = Sp n /£ n belonging to the 
nth step in the energy cascade, one can estimate accelerations by |a| = linin^ooan. 
For high Reynolds number Re 3> 1, or for the situation where effects of the kinematic 
viscosity v can be neglected compared with those of the turbulent viscosity, the 
Navier-Stokes equation is invariant under the scale transformation [ [2U IT5] : 
r — > Ar, u — > \ a / 3 u, t — > \ 1 ~ a / 3 t and (p/p) — > A 2ci//3 (p/p)- The exponent a is an 
arbitrary real quantity which specifies the degree of singularity in the acceleration 
for a < 1.5, i.e., lim^oo a n = lim£ n ^o Sp n /£ n ~ lini£ n ^o ^ 2a/ ' 3 ' )_1 — > oo which can be 
seen with the relation Sp n /Spo = (£ n /£o) 2a ^ 3 - 

The multifractal analysis rests on the assumption that the distribution of the 
exponent a is multifractal, and that the probability P^ n \a)da to find, at a point in 
physical space, an eddy of size £ n having a value of the degree of singularity in the 
range a ~ a + da is given by [ 121 El HI E] 

?W (a) oc [l - (a - a„) 7(^a) 2 l n/(1 ' 9) (2) 



'-mm 



with (Aa) 2 = 2X /[(l — g)ln2]. The range of a is a min < a < a max with a r 
a — Aa and a max = «o + Aa. Here, we assume that the distribution function at 
the nth multifractal depth has the structure P^(a) oc [P^(a)] n . This is consistent 
with the relation [DUE] ?W( a) cx 5 x n *\ a > that is a manifestation of scale inveriance 
and reveals how densely each singularity, labeled by a, fills physical space. Within 
the present model, the multifractal spectrum f(a) is given by [I2HH1SE] 

f(a) = 1 + (1 - q)- 1 log 2 [l - (a - a f / (Aa) 2 } . (3) 

To make the paper self-contained, we put here its brief derivation. The distri- 
bution function (0) is derived by taking an extremum of the generalized entropy, 

lr The velocity fluctuation Su n is denned by Su n — \u(» + £ n ) — where u represents a 

component of the velocity field u. 
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the Renyi entropy [ EDj Sf[P {1) (a)] = (1 - q)' 1 In/ daP {1 \a) q or the Tsallis en- 
tropy [HUH1II3] Sj[pW(ot)] = (1 - (/da pW(a)« - l), under the two con- 
straints, i.e., the normalization of distribution function: / daP^(a) = const, and 
the g-variance being kept constant as a known quantity: <r 2 = (/ daP <yl \a) q {a — 
a ) 2 )/ / daP( 1 \a) q . In spite of the different characteristics of these entropies, the dis- 
tribution function giving their extremum has the common structure (0). The depen- 
dence of the parameters ato, X and q on the intermittency exponent fi is determined, 
self-consistently, with the help of the three independent equations, i.e., the energy 
conservation: (e n ) = e, the definition of the intermittency exponent fi: (e 2 ) = e 2 S~^, 
and the scaling relation: 2 1/(1 — q) = l/a_ — l/a+ with a± satisfying f(a±) = 0. 
The average (• • •) is taken with P^(a). For the region where the value of fi is usually 
observed, i.e., 0.13 < fi < 0.40, the three self-consistent equations are solved to give 
the approximate equations [ 8 : a = 0.9989 + 0.5814//, X = -2.848 x 10" 3 + 1.198/x 
and q = -1.507 + 20.58/i - 97.ll/i 2 + 260.4/i 3 - 365.4// + 208.3//. 

Since there are two mechanisms in turbulent flow to rule its evolution, i.e., the 
one controlled by the kinematic viscosity that takes care thermal fluctuations, and 
the other by the turbulent viscosity that is responsible for intermittent fluctuations 
related to the singularities in acceleration, it may be reasonable to assume that the 
probability K^ n \y n )dy n to find the scaled pressure fluctuations \y n \ = 5p n /5po in the 
range y n ~ y n + dy n has two independent origins: 

M n) (y n )dy n = AP(\y n \)dy n + AA<- n \y n )dy n . (4) 

The singular part Ag^(|y n |) of the PDF stemmed from multifractal distribution of 
the singularities, and the correction part AA^ n \y n ) from thermal dissipation and/or 
measurement error. 3 The former is derived through h^\\y n \)dy n oc P^(a)da with 
the transformation of the variables: \y n \ — 5 2a ^ 3 . The mth moments of the pressure 
fluctuations are given by 

/CO , , 

dy n \y n \ m ^ n \yn) = Vy£> + (1 - 2#) a 2m 5^ (5) 
-co 

with a 3 g = {2/[y/C q {l + y/Q)]} 172 , Cg = l + 2q 2 (1 - q) X In 2 and 

/CO 
dy n \y n \ m A^\y n ). (6) 
-co 

We used the normalization: ((1)) = 1. The quantity 

«nm 2Xm 2 1 r , / / \i / s 

^-KW^"^ [l " log2(l+V ^ ] (7) 

2 The scaling relation is a generalization of the one derived first in [IJSJIJB] to the case where 
the multifractal spectrum has negative values 

3 Needless to say that each term in Q is a multiple of two PDF's, i.e., the PDF for one of the two 
independent origins to realize and the conditional PDF for a value u> n in the range tu n ~ u> n + dw n 
to come out. This is of course in a generalized sense in which the second correction term may 
weaken the first singular contribution. 
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Figure 1. PDF of accelerations plotted on (a) log and (b) linear scale. Compari- 
son between the experimentally measured PDF of fluid particle accelerations by 
Bodenschatz et al. at R\ = 970 (Re = 62 700) and the present theoretical PDF 
(u} n ). Open squares are the experimental data points on the left hand side 
of the PDF, whereas open circles are those on the right hand side. Solid lines 
represent the curves given by the present theory (JHJ) with fi = 0.259 (q = 0.431) 
and n = 16.4. 



is the so-called scaling exponent of velocity structure function, whose expression 
was derived first by the present authors [0E1II1E]- The formula explains quite well 
experimental data [ El El IH El El IB] ■ Note that the formula is independent of the 
length i n , and, therefore, independent of n. 

Let us derive the PDF A^(uj n ) of the accelerations by dividing it into two parts 
with respect to u n , i.e., 



A (n) ( 



uJ.„ 




for \uj n ^ 
for ut < 



< col 



^ , .max 



(8) 



where the scaled variable u n is defined by \u n \ = a n /((a 2 )) 1//2 = cov^ 2 "/ 3 "> 4//2 with 



u n = [2#<^ + (1 



27Q™^)a 4 ] 1//2 , and tu™ ax = cov^ 2 "™/ 3 ?4 / 2 . Assuming that, for 



smaller accelerations \u n \ < u^, the contribution to the PDF comes, mainly, from 
thermal fluctuations related to the kinematic viscosity or from measurement error, 
we take for the PDF A<}(oj n )uj n = Ag l \y n )dy n + AA^ (y n )dy n a Gaussian function, 
i.e., 



A 



(n) 

s ex P 



1 + 



1 



(9) 



with Ag 1 ' ) = 3(1 — 27g n '')/(4u; n y / 27rX| ln<5 n |). On the other hand, we assume that 
the main contribution to Ai< (oo n ) ma y come from the multifractal distribution of 
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singularities related to the turbulent viscosity, i.e., A&(L) n )du n = Ag 1 ' (\y n \)dy n : 



A 



(n) u n 



UJr, 



1 - 



1 - q {3\n\u n /u nfi \) d 



71 



8X1 In <L 



n/(l-q) 



(10) 



with |w ft)0 | = uJl 00 ^-^ 2 . The point was defined by w\ = u n 8^ l z ~M 2 with a \ 
being the smaller solution of C 4 /2 — 2a/3 + 1 — /(a) = 0, at which M n '(uj^) has the 
least n-dependence for large n. Here, A<}(u n ) and A±<(a> n ) were connected at cjt 
under the condition that they should have the same value and the same slope there. 
The specific form of Gaussian function comes out through the connection. 

With the expression (jHJ), we can obtain AA^ n ^(?/ n ), and have the formula to 
evaluate 7^ in the form 



where 



Tin) 



27^ = - Lfc)) / (l + - I$> 

^2(m+l) a t/3_ C4 / 2 ! 

u r> I 



2^27TX\lnS n 

3 ^ mat / 3 
2 A /27rX|hi(i, 



/ ciz z m exp 

JO 



1 + ~/V) 



m— 1 



(3 In |z/4l 



n 



8X1 In 5. 



(11) 

(^ 2 -l)}, (12) 

»/(!-?) 

(13) 



with z n 



= 52(a ma x-at)/3 &nd ^ = ^/^t = ^(ao-atJ/S^ Now> the p DF 

of fluid particle accelerations (JHJ) is determined by the intermittency exponent \x and 
the number n of steps in the energy cascade which gives the eddy size £ n . 

The comparison between the present PDF of accelerations and that measured in 
the experiment [I2T|I21Z] at R\ = 970 is plotted in Fig. [TJon log and linear scale. The 
intermittency exponent \x = 0.259 and the number n = 16.4 of steps in the cascade 
are extracted by the method of least squares with respect to the logarithm of PDF's 
as the best fit of our theoretical formulae (jHJ) to the observed values of the PDF [ 
El 122] by discarding those points whose values are less than ~ 2 x 10~ 6 since they 
scatter largely in the log scale. Substituting the extracted value of fi into the self- 
consistent equations, we have the values of parameters: q = 0.431, an = 1-149 and 
X = 0.307. With these values, other quantities are determined, e.g., Aa = 1.248, 
a + — «o = oiq — a_ = 0.7211 and uj u = 0.557 (a* = 1.005). We see an excellent 
agreement between the measured PDF of accelerations and the analytical formula 
of PDF derived by the present self-consistent multifractal analysis. The value of ui^ 
tells us that the contribution from thermal fluctuation and/or measurement error 
is restricted to smaller values of ui n , i.e., less than about one half of its standard 
deviation. Then, the values of a responsible for the intermittency due to the scale 
invariance turn out to be smaller than a'«l. This is the range within the condition 
a < 1.5 in which the singularity appears in fluid particle accelerations. 

The values of the scaling exponents ( m , given by (J7J), for m = 1 • ■ ■ 20 are listed 
in Tabled for future convenience in comparison with experiments or other theories. 
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Table 1. The values of the scaling exponents derived by the present theory for 
the acceleration is also listed with the corresponding values for q, ao and X. 



m 


Cm 


m 


Cm 


m 


Cm 


m 


Cm 


1 


0.3661 


6 


1.741 


11 


2.599 


16 


3.175 


2 


0.6989 


7 


1.945 


12 


2.731 


17 


3.268 


3 


1.000 


8 


2.130 


13 


2.854 


18 


3.356 


4 


1.272 


9 


2.299 


14 


2.969 


19 


3.438 


5 


1.519 


10 


2.455 


15 


3.075 


20 


3.515 



Note that /x is related to C6 by the relation n = 2 — ( e within the present analysis. 
The flatness = ((a^))/((a^)) 2 = ((a>£)) of the PDF of accelerations has the value 
F^ n ' = 57.8 which is compatible with the value of the flatness ~ 60 reported in 
[ |2H 122] as it should be. It is quite attractive to see that the distance r = £ n 
corresponding to the extracted value n = 16.4 reduces to r = 0.821xxm (r/rj = 
0.0456), which is close to the value of the spatial resolution 0.5/xm (1/40 of the 
Kolmogorov distance) of the measurement in [ [^J |22] • 

We derived, in this letter, the formula of the PDF of accelerations (JHJ) within 
the approach of the multifractal analysis constructed by the present authors, and 
analyzed successfully the beautiful experiment conducted by Bodenschatz and co- 
workers [1^11221 m the Lagrangian frame. We expect that other data for the same 
system, in addition to the PDF of accelerations, will be provided in the near future, 
such as the intermittency exponent /x, the scaling exponents Cm; the PDF's of the ve- 
locity fluctuations or of the velocity derivatives. Then, we can cross-check the validity 
of the present analysis based on the formula (jHJ)- 4 Note that the formula for the PDF 
of accelerations is different from the one for the PDF's of velocity fluctuations or of 
velocity derivatives. The empirical PDF, A emp (c<j) = C exp {— u 2 / [(1 + \ojf3/a\ y ) a 2 }} 
with f3 = 0.539, 7 = 1.588, a = 0.508 and C = 0.786 proposed in [EQE2I for the 
data at R\ = 970, gives a line very close to the one provided by the present PDF 
A( n )(u;) of (JSJ) for the region \cu\ < 30 where the data of PDF exist. They deviates, 
however, for \cu\ > 30, i.e., A cmp (cu) < M n \u). The extraction of the PDF of accel- 
erations out of the DNS data obtained by Gotoh et al. [ |2Zj is one of the attractive 
investigations in order to check the validity of the present formula (jSJ) derived by 
the unified approach providing various PDF's [ El HI M M j since the accuracy of 
their DNS data for PDF's is very high up to the order of 10 -10 [I28j. 

The authors would like to thank Prof. C. Tsallis for his fruitful comments with 
encouragement. The authors are grateful to Prof. E. Bodenschatz for his kindness 

4 We have checked, with the help of the DNS data reported by Gotoh et al. in [H7|, the accuracy 
of the value of the intermittency exponent extracted out of the measured PDF's. With the formula 
0, we determined in [|Hj the value /x = 0.240 for the longitudinal velocity fluctuations by fitting, 
with the method of least squares, the ten data of the scaling exponents £ m (m = 1, 2, • • • , 10) at 
Ra = 381. On the other hand, we have extracted the value fi = 0.237 by the method of least 
squares as the best fit of the theoretical formula for the PDF of velocity derivatives, derived in [ 
EI by the multifractal analysis, to the corresponding PDF data. 
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to show his data prior to publication. 
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